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REAL TIME DIGITAL PROPULSION SYSTEM SIMULATION 
FOR MANNED FLIGHT SIMULATORS 
by James R. Mihaluew und Clint E. Hart 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 

ABSTIUCT 

In Uie past propulsion system simulations used in High', simulators 
have been extremely simple. ThisTvsults in a loss of simi lation realism, 
eliminates significant engine and aircralt interaction dynamics and pre- 
vents generation of important internal engine parameters. Reasonably 
detailed simulations will be necessary to jjermit system evaluations in 
a simulated flight environment. A real time digital simulation of a STOL 
propulsion system was devolo|)ed which generates significant dynamics 
aixi internal variables needed to evaluate system iierformance and air- 
craft interactions using manned flight simulators. The simulation ran 
at a real-to-execution time ratio of 8.H. The model was used in a piloted 
NASA flight simulator program to evaluate the simulation technique and 
the propulsion system digital control. The simulation is described ?nd 
results shown. Limited results of the flight simulation program are 
also presented. 


INTRODUCTION 

Cost considerations are becoming more importart in planning ex- 
perimental and flight programs. Simulation, with its inherent analytical 
flexibility, can be used effectively to develop aircraft and propulsion 
systems to a higher degree before hardware is built and tested. It also 
provides a safe means of studying failure modes and effects. 
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Studies usin^ niutmod simulators to evaluate fli(;ht stability 

and enj;ine-out ix'rformanee have iH'en made for STOL ty|X‘ aircraft 
(1,2). Propulsion system models used in these simulation studies are 
extremely simple providing only a tiirust sit^nal. No concern was t;iven 
to how Uie thrust was derived or what mifjht bt' hap|X'nin^ to Uie entjine 
producing tJuit tlirust. For example, the engine simulation used in (1) 
and (2) was derived from a series of thrust transients from an early 
oTOL engine hybrid com|XJter model. The res|X)nse was matcluKl to a 
transfer function which gave a typical resixmse. 

Fngine models were also kept simple Ix'cause of the limited com- 
putation time available. Using them' imxleL'} resulted in a loss of simu- 
lation realism and eliminated significant engine and aircraft interaction 
dynamics. Also, ini|X)rtant internal engine parai.i^ters were unavailable 
for analysis. To overcome Liese deficiencies, reasonably detailed real 
time propulsion system simulations must Ix' develojx'd. Tliesc will pro- 
vide the capability to evaluate propulsion systems and Uieir interaction 
with aircraft controls on a real time basis in a simulated flight environ- 
ment. 

The QCSKE (Quiet, Clean Short -liaul Exix'rimental Engine) Program 
was initattHl by NASA to develop and demonstrate iiropulsion system tech- 
nology for an advanced commercial STOL aircr;ift. One of the s|H»cific 
technical objectives was to provide technology for digital electronic con- 
trol of future commercial engines. An element of this technology devel- 
opment was to evaluate tlie digital control in a simulated flight environ- 
ment. 

To accomplish tJie QCSEE Program requirement an effort was 
initated whose overall pur|X)se was to evaluate the QCSEB] UTW (Under- 
tlie-Wing) digital control system over a range of conditions encountered 
in typical air|x)rt oix'rations. The goal of Uie simulation effort was to 
derive a real time digital propulsion simulation which could be integrated 
into a multi-engine aircraft simulation. 

Tills iKiper summarizes Uie accomplishments to date. It is divided 
into two parts. F'irst, the simulation model and tecliniques are described 
and results presented. Second, the application of tlie propulsion system 
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simulation to an aircraft simulation is presented. Limited results of 
the flight simulation test program are presented. 

REAL TIME DIGITAL SIMULATION REQUIREMENTS 

The nature of the digital comixitipr, process to ijerform real time 
digital simulation is one of efficient mathematical modeling to attain 
the desired detail \uth minimum computation time and maximum numer* 
leal stability within u s|)ecified time |)eriod. In this study we wish to 
model tile propulsion system to a degree such that its internal engine 
and control system parameters will be available for analysis. For the 
level of steady state and dynamic complexity required to meet program 
objectives, steady state accuracy does not have to be compromised over 
detailed models. However, high frequency content may be reduced sig- 
nificantly. The initial improvement to Uie simulation comes from effi- 
cient programming to attain minimum calculation time. More rapid 
improvement occurs by maximizing numerical stability to i>ermit long 
time steps. 

The comixiter cannot meet real time requirements simply by re- 
{x?atir.g Uie calculations as rapidly as |X)ssible. The time net>ded to 
make the calculations bears no relation to the time step used to change 
the time variable. A means to schedule and account for computation 
time is essential to assure that Uie computations can be accomplished 
in the allotted time. Time available for computation of Uie propulsion 
system must be shared since it is only a part of a larger simulation. 
Real time, then, in the context of overall simulation implies that the 
propulsion system will have to be calculated faster than real time to 
be effective. 


PROPULSION SYSTEM SIMULATION 

In Uiis part of the pajier, the prc.pulsion system is described and 
the development of the simulation is presented. 
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QCSKE Ex|)oriniontul Propulsion System 

Engine . - The UTW engine is shown in figure 1. The engine uses 
an FlOl core gas generator with a high bypass fan duct. The system 
features a higti Macli inlet, a variable pitch fan, a variable geometry 
fan duct oxliaust nozzle and a digital electronic control system com- 
bined with a hydromechaiiical fuel control. It is designed to provide 
'^7395N (17400 lb) of iiistalled thrust at takeoff on a 305.6 K (90 F) day. 

Tne fan is a low pressure ratio, low tip s|X'ed configuration witii 
variable; pitch blades and is driven by a low pressure turbine tlirough 
reduction gears. The fan is capable of blade pitch changes from for- 
ward to reverse thrust. The fan variable pitch actuation and control 
is designed to move the blades from forward to reverse {xisition in less 
than one second. 

The fan exiiaust nozzle is a hydraulically actuated variable area 
design. It is capable of area change from takeoff to cruise as well as 
oix'ning to a flare |x>sition to form an inlet in the reverse thrust mode. 

Control. - The control system manipulates four variables to achieve 
rapid Uirust res|x>nse and noise suppression. Control of engine pressure 
ratio, fan s|)eed and inlet Mach number is accomplished by manipulating 
fuel flow, fan blade pitch and exiiaust nozzle area. Variable stator vanes 
are scheduled by core speed to attain optimum stall margin. 

The structuie of Uie control system is shown in figure 2. The digi- 
tal control is the heart of the system and controls the output variables 
in resixmse to commands from the aircraft. It generates all control 
laws and logic and most of the limiting functions as well as ixiwer man- 
agemeiil, condition monitoring and failure corrective action and indica- 
tion. The hydromechanical control provides an electrohydraulic servo 
fuel valve which is used by Uie digital conU'ol for primary fuel control. 

It also provides backup fuel control through a core sjx'ed controller, 
acceleration and deceleration limits and primary control of tlie core 
compressor stators. The hydromechanical throttle is used as a mech- 
anical enable or power limiter. 
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Analytical Model 

Tlie aiuilytical model is derived Irom the real QCSEK UTW pro- 
pulsion system. It represents mathematically the steady slate and 
dynamic relations between the eii(;ine com|x>nent representation and 
the control coni|X)nent representations. Engine dynamics are based 
on the dynamic form of the conservation equations and engine transient 
exiJerience. Steady state engine |X‘rformance is based on com|X)nent 
representations derived from cycle model data and steady state forms 
of the conservation equations. The variable pitch fan and its ability to 
reverse duct flow imixise unique engine modeling problems. 

Engi ne repre senta tion . - A major initial assumption is to eliminate 
high frequency elements since Uiey are not necessary to meet simulation 
objectives. As a result only rotor dynamics and compressor and turbine 
heat capacitances will be retained as true states due to their low fre- 
quency content. Also as a result of this assumption, certain iterative 
variables will be required. All other states in the engine are neglected 
and Uie component ix'rformance maps are manipulated functionally to 
accomodate these assumptions. 

The form of the engine model and information flow follows the sche- 
matic shown in figure 3. All major engine components for the bypass 
duct and core are represented including the inlet throat and duct |x»rfor- 
mance. Tliese are im()ortant since the QCSEE UTW engine oijeratet 
witii a high Mach inlet. 

Fan tip ixrformance is represented using maps for pressure and 
teni|)erature ratio as functions of corrected srx?ed and weight flow. In 
this case tlie variable pitch angle is added as a parameter. The fan hub 
|)erformance is assumed correlatable to the fan tip jX'rformar.ce. Coie 
inlet duci loss is also included. Fan inlet and discharge pressure are 
chosen as iterative v.ariables. Fan duct pressure loss is derived from 
pressure loss relations as functions of duct airflow. The fan noz.’.le is 
represented by the nozzle flow equation for unchoked flow. Nozzle area 
is variable and constant discharge and velocity coefficients are assumed. 

The compressor is represented by maps for pressure and tempera- 
ture ratio. It is assumed that bleed flows are constant iiercentage of 


6 


Inlet flow ;md that Ute compi'cssur o|)crate8 with tiie variable stator 
vanes on schedule at all times. Compressor heat capacitiince dynamics 
are derived from experimental correlations. 

The combustor is represented by relations which include pressure 
drop and heat rise. Combustor dynamics are nei^lected. Combustor 
dischari;e pressure is chosen as a required iteration variable. 

The h4;h and low pressure turbine are represented by flow and 
enthalpy drop maps as functions of pressure ratio and corrected s|X'ed. 
Heat capacitances are derived Irom experimental ex|H'rience. 

The core nozzle is represented as the fan nozzle except that area is 
constant. The nozzle inlet pressure is a required iteration variable. 

Control representatio n. - Since the control evaluation is a prime 
objective of the profjram a detailed control representation is essential 
es|)ecially for tlie control laws and switching logic. 

The controls model includes representations for the digital a.id 
hydromechanical controls conqionents. The digital portion contains 
detailed representations of Uie fan sjx'ed, inlet Mach numlx'r and en- 
gine pressure ratio (jiower) controls. The hydromechanical control 
includes Uie core s|x>ed, acceleiation and deceleration controls. 

The schematic of the control system is shown in figure 4. This 
diagram will not be discussed in detail since it is s|)ecific to Uie QCSEE 
engine. It is Included to illustrate the extent of the control representat- 
tion ii; Uie simulation. For further details see (3) and (4). 

Tile resulting luialytical model is sixteenth order and includes four 
engine states, four iterative loops, four control sensor states and four 
control actuator states. These do not include the digital controller 
suites. 

To meet Uie real time requirement, meUiods were required to gen- 
erate the functional form of the coni|Xjnent representations and to attain 
convergence in Uie iterative variables. These techniques are described 
next. 
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Futu’ t toil IU*pri‘ 80 iit at ion 

Kxaminution of llu' fuiu lion roproHonlalions roquiriHi in ihiH niodc‘1 
8how that 8onu' art* niultivariai)U> fuiu tionn ol two or tliroe variatiU'H. 
To inipU'iiuMit thoHo reprosontationH in tiu' Hiinula'ion, two (xilntH niuHt 
Ih' ronHidert'd. First, tho fuiu'tioiu. niuHl rcprt'm'nt tIu' |H‘rforniaiu'(> 
of tlu' en(;iiu' coni|H)iu'nlH over a wide rani’e of o|X'ratiiai. SeeoiKi, f'.ie 
eoni|xitationaI method develo|XHl must Ih> accurate and fuHt to meet the 
real time Hiieulation requirement. 

Functi on si mpi ific atmn. - A detailed QCSKF cycle model proi;ram 
wau us(h 1 to obtain data for tlie comixment function repreHentatioiiH. 
Since the fan ujX'rateH over a wider ran^e of conditions than the cum* 
pressor or turbines, various combinations of corrected fan siieiKl, 
corrected weinlit flow and pitcli an^le were selected as inixits to tlie 
pro(i;ram. The out|xits of the proi;ram consisted ol the steady state 
values of most of t!ie eii(;int internal variables. 

Fro.n this data the iH'rformance of eacli comjxment was plotted 
according to tJie functional relations required in tiie nuxiel. In this 
way tlie o|X'ratin(; ranije of eacli com|xment in relation to Uie fan o|X'r- 
atinp ran^e mas determined. A numlHM* of simplUications became 
evioent. The most si|;nificant of tiiese was that the cure flow patli was 
choked over tJie entire fan o|x>ratini; ran^e. The compressor and hi|jh 
{n'essure turbine mails could Uius 1 h' reductni to functions of a siiii;le 
vai iable. 

Function t;enerator s. - Many iteneralized function generation 
routines for single and multivariable functions have been develo|)ed. 
Tiicse are based on table search and inter|X)lation methods. Function 
generation routines descriluxl in (5l are convenient to use and calcula- 
tionally elficient. Tliese routines were usihI to program the function 
generation requirements of the simulation. Tests were made of the 
calculation time required for various tasks in the simulation. It was 
found that approximately 40 ix'rcent of the total calculation time was 
consumed in function generation. Tins indicated that calculation time 
could Ih' ilei reased significantly if faster function generation methods 
could be found. 


Curve fitting . - AnoUier simple, direct means to accomplish func- 
tion representations is curve fittini;. Curves representlnn single or 
multivariable functions can Ix' matched by equations containing poly- 
nomial. geometric or other analytical functions. Comixiter programs 
are available to carry out the details of various curve fitting techniques; 
however, much trial and error is Involved. Higher order equations are 
used to Improve accuracy but comixitation time increases. Higher accu- 
racy can also l)o attained by dividing tlie curves into segments and match- 
ing each segment with lower order equations. The curve fitting techni- 
que was eventually adopted ior Uiis model development. 

To illustrate this function generation method, curve fitting techni- 
ques are used here to develop tJie fim corrected flow map. The fan 
corrected flow map is a function of three variables: pressure ratio, 
corrected sjx'ed and blade pitch angle. It was assumed that the fan map 
could be represented by the functional relation shown in figure 5. Th<s 
representation consists of a basic fan map which is a two-variable func- 
tion based on data for zero-degree blade pitch angle. Modifying func- 
tions are applied to one input and the output of Uie basic fan map. 

Using data obtained from the QCSEE cycle model program Uie func- 
tions were plotted to determine Uieir functional form. The S|)eed lines 
of the basic fan map could be represented by hyix'rbolas. The coeffi- 
cients of the hyix'rbolas were functions of corrected s|X'ed calculated 
by third order polynomial equations. The pressure ratio function was 
found to be sensitive to both corrected s|x>ed and pitch angle. It was 
represented by a function of two variables derived from a multiple 
regression formula. Tlie flow function was sensitive only to pitch angle. 

It was divided into four segments represented by two linear and two 
second order polynomial equations. 

A comparison is shown in figure 6 for the basic fan map cycle model 

data and the generated function data. The comparison is excellent. The 

map generated at a constant corrected speed of 95 percent and various 

fiui pitch angles is shown in figure 7. The accuracy is reasonably good 

* 2 2 * 

for negative pitch angles along the constant 1.613 m (2500 in ) nozzle 
area line where most of the oix>ration occurs. More accuracy could be 
attamed but only at the cost of more complexity. 
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Tht* imuiuKl of function t^cMicration liantHi on curve tttting lochniqucH 
proved to Ih‘ cnlculationally KiMter than the iieneralized metluid for all 
functloiu*. The lai'i^eHt rt*duction in calculation time occuru for multi* 
variable funettonn. For the three variable fan map the calculation lime 
|H>r |)oint wan 85 pn UHlni; the curve fiitini: method. The i^eneraltzinl 
function iieneratlon metiuMi r(*qulrtHl 750 pH. 

All com|xmenl repreHentation were correl.iied uHin»i this icilinique 
It eliminateH all generalized function geiM'ralion and Hubroutine callH. 
liH main disadvantage at this |M)int iH its lack of gt'iieral application and 
associated time consuming optimi/.ation using standard curve fitting 
routines. 


Convergence MelluKls 

At this |K)inl in the nuKlel development the functional relations have 
iH'en determined ami the states esiablisluHl. In addition to the states 
which require integration, the four iterative Uxips require a metluKl to 
converge to a solution. 

Iterat ion. - The initial convergence technique used was a simple 
iieration procedure. The equation formulation and solution order re- 
sulttHi in one double nested loop and two single liH>p iterations TIu* 
basic pimblem with iteration was that tlie lime to converge was tran- 
sient depemient; that is, cycle time varied pro(X)rtionally to the rapidity 
of Uie transient. This is to Ih» extx'cled since larger nonsteady mis- 
matclu's occur during transients than near steady slate and more itera- 
tions are requinxi to converge the difference. Simulation frame time 
must be based on the maximum calculation time in order to avoid missed 
intervals. Kfficlency of this metluxl is lud goixi and was considertni un- 
acceptable in this nuxiel development. 

Integration . - Litegration is anotiier metluxi which may Ih' used to 
provide approximate convergence in iterative loops. In tliis application 
ihe iteraiion differences which occur are inte;:rated. A sketcli of the 
process is shown in figure H. The procixiure is similar to using high 
gain integrators in analog comimiation to prevent algebraic loo|>s. The 
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previous "guessed*' value, Is used to calculate the new required 

steady state value, P. The difference between Uiese two values is inte- 
grated to generate an ufxiated "guessed" value. The process shares 
the same integration algorithm used for Uie model states. 

Integration algorlUim evaluation. - To select a suitable integration 
method an evaluation of several one step direct algoritlims was made 
using Uie model equations. Tlie fourth order Hunge-Kutta method was 
used as a standard measure of accuracy and stability. On the otlier 
end of Uie six'ctrum, simple Kuler integration was also evaluated. No 
error analysis was made. Instead an empirical approach was used to 
determine ix'rformance based on a standard model transient. 

Since four derivative evaluations are required |)er time step for 
the Hunge-Kutta, accuracy and stability are high but calculation time 
is also high even thougli laifjtr time steps are ixissible. On the other 
hand tlie Euler metliod Itad relatively large errors and stability was 
|X)or for Uie same time step range. 

Two second order methods, Uie improved Euler (slo|)e averaging) 
and Uie modified Euler (|)oint averaging) metluxis were also evaluated 
and found to be comparable. The modified Euler had a slight advantage 
in programming efficiency. It also provided a compromise in Uiat it 
coni|xited in half the time of the Hunge-Kutta but provided only slightly 
less stability. The algoriUim proved to be stable over a wide range of 
transients. Because of these advantages the modified Euler algorithm 
was selected for use in the real time simulation. 

Simulation Model 

Using the functional correlations derived )or the engine comixments, 
the state equations, the digital control algorithm and tlie integration con- 
vergence technique, a simulation model was assembled and programmed. 
Tlie model resulted in a set of sixteen simultaneous first order nonlinear 
differential equations. These are solved as an initial value problem. 

Every effort was made in programming to eliminate unnecessary 
calculations. Divide compuiations were minimized and exponentiation 
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avoided. KxIonHivp um* of branching; on condilion wa8 also used. 

The result in |4 rom|Hiler program (low diagram is shown in fi(;ure 9 
The model programming and development was aecomplished on a 
Univae 1100-40 computer using Fortran V On this computer the 
(light segment of the model consumes 1.9 ms cominitation time |x>r 
time increment. Cycle time was not transient sensitive On the 
flight simulator computer, a Xerox Sigma 8. cycle time was 5.7 ms 
or three times slower. On both com|XJters the stmulatton was found 
to be stable and accurate for integration lime steps up to 50 ms. 

Thus the real time lo-execution time ratio was 26 for the Univac and 
8.8 for the Xerox. 


Simulation Ifesults 

Simulation verification. • No ex|M<rimental transient data on the 
real engine and control system arc available to verify the model re- 
sults inde|x*ndently. Since the detailed cycle detk of the engine is the 
only source of iix1e|X‘ndent steady state data, a com|)ari8on is made 
with It. 

Figure 10 shows a comparison of cycle deck values and simulation 
model data for a numlx>r of selected parameters at o|)erating |X)lnts 
which were of primary tm|xirtance (or the control evaluation. These 
weie; 100 |X'rcent takeoff jxiwer, 62.5 |X‘rcent approach |X)wer and 
100 |x>rcent reverse fxjwer at sea level standard day conditions. The 
comparison data were generated by setting simulation control values 
of |X*rcent corrected fan speed, fan pitch angle and fan exhaust nozzle 
area into the cycle deck. The cycle dec k data, model data and |x*rcent 
error are shown Parameters are defined according to the engine nnxlel 
schematic in figure 3. 

The simulation data shows excellent agreement with the cycle deck 
for Uie two forward cases Errors are generally less than one iJercent 
for the parameters listed. The reverse case is the least comparable 
with errors in the core path slightly over 5 |)ercen». The mixiel steady 
state accuracy is very acceptable. 
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Transient performance . - A number of transients were run which 
normally can be ex|)ected in a manned flight simulation program. These 
included normal as well as failed forward and reverse transients. A 
limited numiH>r are presented here to indicate the model's flexibility 
aiKl to show the detail and extent to which engine variables are available. 

Normal forward transient: One of Uie specific QCSEE program 
objectives was to develop a control which would insure a transient re- 
s|X)nse of less than one second from 62.5 to 95 percent net thrust. 

This is accomplished primarily by maintaining a higli fan speed witli 
manipulation of Uie fan pitch. Fuel is manipulated to maintain a sched- 
uled engine pressure ratio while the exhaust nozzle is o|)en to a high 
area setting. 

Figure 11 shows a typical approach waveoff transient from approach 
|X)wer at 62.5 iiercent to go-around |X)wer at 100 jiercent. Shown are a 
number of engine and control parameters as a function of time. Practi- 
cally ' M • ngine or control iiarameter can be extracted from the model 
for i;!splay. These shown are only representative. 

As shown in the figure th.' required thrust is achieved in about 0.6 
oeconds. As |X)wer is stepiied tlie power control selects the turbine inlet 
tem|x.‘rature limit loop and then the acceleration limit and subsequently 
switches back to Uie engine pressure ratio control loop. The fan exhaust 
nozzle moves rapidly to an area slightly less Uian the takeoff value and 
Uien moves to control inlet Mach number for noise suppression. Like- 
wise, Uie fan pitch is rapidly changed to maintain fan s)x?ed and, as 
shown in the fan sjK'ed transient, is a contributing factor in producing 
Uie rapid Uirust increase. 

Normal reverse transient: AnoUier of the specific program objec- 
tives of QCSEE is to provide for a fast Uirust reversal capability. A key 
factor here is in Uie design and evaluation of Uie control system logic. 
The thrust reversal must be achieved in less Uian 1.5 seconds while 
maintaining safe engine oiieration. The variable pitch is used to reverse 
direction of the fan duct airflow. In Uiis model, reversal occurs by 
rotating the blades through stall. 

Figure 12 shows Uie transient res|xmse of a simulated takeoff abort- 
to-reverse sequence. As shown in the transients, U.rust decays rapidly 



tj a roverHo U'vi'l. riu> small Ihriist iiuToaso a( reverse initiation is 
due to a flow inerease priKlueed by the Ian piteli o|X'ntii|; iH'lore the fan 
S|K>ed res|)t)nds. The flat sef^ment at zero thrust is an effect simulating 
flow reattaehment as the fan blades i^o lhrou(;h stall. The same effects 
occur in com mi; out of reverse. A steady state reverse thrust level Is 
approached in alxuit l.d seconds. 

C omp uter failure transient: Figure 13 indicates the transient whic'i 
occurs when a simulated computer failure is im|X)sed. The failure is 
simulatt‘d by settini; the dii;ital computer servovalve torque motor cur- 
rents to zero. In this condition the nozzle i;(Hf>8 to a wide o|M'Ii |X)sition, 
the fan blade pitch fails fixed due to a faiUfixiKl servovalve and fuel 
control IS assumed by the hydromechanical core S|H*ed control. The 
core SjX'ed scliedule is set to 10 jw'rcent alxive that which would occur 
durinp: normal steady state o|X'ration. 

Die initial |X)wer is at 02.5 ix'rcent. This is followed by a step 
increase in |>ower to 00 |H‘rcent. The com|Kiter failure occurs at 5 sec- 
onds. As shown, net thrust l.icreases as well as fan and core S|x*ed. 
Control IS still maintaim'd through the |X)wer setting; input. 

Other transients: Other normal and failed transients were *-un 
includinf; takeoff rolls to altitude, constant throttle climb to altitude, 
inadvertent icv?rse on takt'off and apfiroach, and various servovalve 
failure modes. It is important to inv''sti^;ate every conceivable condi- 
tion that may be encountered in a simulated flii;ht evaluation primarily 
to insure safe flif;ht simulator o|X*ration and to determine model ojx'ra- 
tional limits. 


FLIGHT SIMULATOH APPLICATION 

This .sectmn describes the application of Uio propulsion simulation 
to a ilit;ht simulator pro(;ram used to evaluate the simulation techniques 
aixl the propulsion system control. The simulation program used tJie 
NASA Ames Flij;ht Simulator for Advanced Aircraft (FSAA) facility. 

The QCSFF UTW simulation develo|X‘d previously was interfaced with 
the Fxternally Blown Flap (KBFj airframe (O). Tlie model integration. 



14 


( 


aiwl t'vaUiadon program arc (U‘HcriU>d ami typical prclimi* 
rury i chuUh arc prcHcntiHl. 

Flight Simiilatiir Description 

The FSAA is a i;(MMM'al |Hir|M>s<‘ aircralt simulator designed for 
trans|)ort aircraft research. It uses a two>man simulator cockpit 
o|H‘ratinp in six dep;rees of motion with lateral motion u.) tii 100 feet. 

It IS supiHM'tiHl l>v a central computing facility and a terrain subsystem 
for visual efl«>cts. A briel di>scription of tlu* facility is (iiven liere. 

For a more detaileil dt'scnption ami user’s ^uide see 17). A com|Hisite 
photo ol die system is shown m fi^urt' 14. 

Computer . - I'he real time simulation facility includes a digital 
com|Hiter will) links to all |M't*i|)iieral eipiipmeni relevant to real time 
simulations. The hi'arl is llu' digital computer. It is a Xerox Si^ma 8 
with a 128 K word core usin^ 82 hits |M'r word. It has a cycle time of 
2.5 /is. 

Sofiwari'. • Ift'al lime simulation of virtually any aircraft nuKiel 
can Ih' implenuMiled by way of a subroutine system. Thest' functions 
include all common kinematic and at'iodynamic relations as well as 
environmental and motion relations. I'lu'.se subroutines are partitioned 
with res|H'cl to frt'quency content into time loops for multirate com- 
putation. 

riie engine simulation is a usi'r supplit'd subroutine. Program 
inputs include throttle |H)silions. failure switches, aircraft states and 
atmosplu'ric pro|H'rlu's. Program outputs include the aircraft axis 
com|xinents of ll'.e applied loripie and forct'S acting on the aircraft due 
to tlu' pro|Hilsion systt'm and paranu ters to drive enft;ine instruments 
in the cab. 

Visua l tl ispla y . - Visual display includes terrain nuHiels which 
have several runway ci>mplexes including a S TOP strip. Ciia|)hics are 
generated from the terram models usin^r a movinj; base camera com- 
manded by the di^;ilal ct>mpuler. 
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Aircraft and Control System Description 

Aircraft . • The aircraft simulation used in this study was a high 
wing four engine externally blown flap (EBF) civil STOL trans|X)rt 
simulation develo|)ed for previous programs. Tlie aircraft was designed 
for use witJi high bypass turbofan engines such as tlie QCSEE UlTV en- 
gine. A complete description including Uie aerodynamic model is given 
in (1) and (6). 

Control system . - The aircraft control system used is described 
in (1) and (6). In addition to the basic aircraft control system, engine 
failure or thrust loss com|)ensation controls were incor|)orated. These 
include a programmed thrust and roll trim coni|)ensation system, a 
thrust command regulated thrust controller and a flight path stabilization 
controller. Characteristics of these control concepts are described in (2). 

Model Integration 

The propulsion system mod<?l described in the first part of this 
paiier was integrated into the EE F STOL airframe simulation. 

Looping structu re. - Normally an aircraft simulation on the FSAA 
is partitioned into three time loops wi!h a 1:2:4 frame time ratio. T!ie 
engine subroutines are also normally called from the slow or low fre- 
quency loop. Frame times for this loop are usually greater than 60 ms 
so tliat the QCSEE models were precluded from lieing called in that loop. 
Additionally, the engine model required more computational time than 
the frame time of the fast loop would {lermit. I herefoi t tlie QCSEE en- 
gine was called from the medium 8|X'ed loop. To save additional com- 
putational time all fast loop subroutines were also called from the medi- 
um s|X'ed loop along with the engine. In addition Uie slow loop was up- 
dated only one-third as often as the medium loop. The resulting loop 
structure is shown in figure 15. Using tliis three-to-one two loop 
structure, enough computation time was available to provide two com- 
plete p..'opulsion system simulations. The entire simulation ran in real 
time on the Sigma 8 with a basic frame lime of 30 ms with the simulation 
IK'rirxi at 150 ms maximuPi. 
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M uiticnt;inc simulation . - Only two coniplcte engine simulations 
were included in the aircraft simulation. Since the aircraft required 
four engines a scheme was devised whereby one simulation was used 
for Uie normal or unf ailed engine and tlie otiier was used for either 
normal or failed oi)eration. A lof;1c subroutine, shown in figure 16, 
etiabled selecting Die engine |x>sition at which tiie programmed failure 
was to oc'cur and placed the failure model at that |X)sition and assigned 
Uie pro|X'r till ottle to it. Control panel discrete signals were used to 
initiate failures at any time during Uie flight. The unfailed engine 
model was assigned to the remaining |x>sitions with tlieir tlirottles 
ganged togeUier for eltlier forward or reverse tiirust o|X'ration. 

In addition a simple engine model, as described and used in (1) and 
(2), was programmed to provide a rapid, simple metliod of trimming 
tile aircraft prior to a flight since the complex models had no fast ini- 
tialization caiiability. The simple nuxlel was matched in steady state 
tlirust, ram drag and engine S|)eed to data from the complex model. 

The simple engine could interface witli Uie EOF airframe in Uie same 
manner as the complex engine models. 

Test Plan 

The simulation was exercised in a test program in which environ- 
ment, aircraft control configuration, propulsion system mode and flight 
mode were vai’ied. 

Environmental variables included wind velocity and direction and 
turbulence magnitude. Aircraft control configurations included un- 
coni|)ensated, thrust command and flight paUi stabilization versions of 
the basic aircraft control. Propulsion system modes included normal, 
abort and failure modes. Flight modes included takeoff, cruise, ap- 
proach, landing and reverse wiUi and without braking action. Failures 
were prescreened and limited to those types that caused pilot concern 
if they occurred in flight. 

In addition to Uie broad objectives stated earlier, a number of 
six*cific objectives were established in Uie test program. These were 
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roluUui (^iMUM'ally to and control iXM lormuncc and pilot reaction. 

Ah rt'(*ardH iH'rformuncc thcHc included lurlunt' life conHuinption. IhruHt 
linearity, Ian blade pitch imHiclalion, reverne control Ui^ic and effect 
of ^UHlH amt turbulence on irlet Mach nunilH'r control. Pilot reaction 
iteniH incliKle IhruHt reH|H)tiHe. en^iiu* aiul control failure reH|ionHe and 
automatic control Ht'quencin^. 

Preiiniinary Kesults 

Data was recorded in (he form of time hiHtoru'H tor ov«‘r 100 tli|;htH 
with various NASA ex|M*ri.nental test pilots. Data recorded include en* 
(;ine |H'rformance variables, control system parameters and aircraft 
status. Output was in (he lorm of multipU'xi'd (two variables |nm’ channel) 
strip charts for selected \ariables on both engine nuHlels. In addition, 
the complete data set of engine, control, environmental and Ciintrol input 
variables were recorderl on ta|H' for each time interval. Also, a print- 
out was obtaiiu'd of seU'cted variables for initial and final conditions, 
minimum-maximum values and statistical variations including; mean and 
staniiard deviations. 

Simulatio n oiHMation al limits. - Ifesults of the simulation ti'chnique 
evaluation are based primarily tm observable restrictions to nuHlel 
o|H'ration during’, the integration and flight test phases. 

Although the propulsion system nuHiel was exercised over (lie ex- 
|H'c(<h1 flndit o|H'ra'.ion ran^e, o|H'rational conditions and engine manipu- 
lations were encountered in the flight program which were not pre- 
examiiu'd. riiese, in general, led to o|H'rational limits which were im- 
post'd on (he model. 

Power range was limited at the law end to lU) iH'rcent in forward 
and 40 jiercenl in reverse. These limits w*'re needed to avoid run time 
errors caused by a nuKiel instability related to the lead-lag representa- 
tion used in the turbine heat soak model. The instalnlity was initiated 
by large rapid throttle decreases which usually occui rtnl after aircraft 
tcuchdown. The only restrictmn to the flight program was tlie require- 
ment to start at a higher thrust level tiuui that of ground idle at takeoff. 


Another bmitation ini]X)sed was the restriction to the flight Mach 
number to a maximum of 0. 4. The thrust parameter used in the QCSEE 
|X)wer control is engine pressure ratio. It is insensitive to fligh: Mach 
number below about 0.4. Above that, however, nonlinear corrections 
to the thrust management schedule must be made to maintain a flat rated 
thrust rating and linear |X)wer -to 'throttle demand. This omission in the 
existing QCSEE control led to tlie restriction on flight Mach number. 

Within the limitations imposed there were no run time errors or 
missed intervals due to the engine simulation during the fligtU simulation 
program. 

Failure effects evalu ation. - A preliminary evaluation of all failure 
modes programmed into the engine simulation was made to determine 
which failures produced effects of consequence to the aircraft. The pur- 
pose was primarily one of safety but also to advise tlie pilots of possible 
effects cn aircraft stability so tliat comiiensatory changes to the aircraft 
control could be anticipated. Failures simulated included: engine out, 
digital computer failure, inadvertent reverse; nozzle, fan pitch and fuel 
servovalve failures in open, closed and fixed {xisitions; and nozzle, fan 
pitch and fuel valve |X)sition and engine pressure ratio sensor loss. 

All failures except those involving the nozzle required significant 
pilot action. The engine out, inadvertent reverse, failed close fan pitch 
and failed closed fuel servovalve failures (in general, failures resulting 
in thrust loss) rc-quired roll trim com|x?nsation. This was provided 
tiirough a conmanded differential ixisition of outboard flap in the lateral 
stability and augmentation system. Normally, this trim is actuated by 
the pilot. However, if a digital control contains failure detection and 
indication logic, an important aspect of digitial control capability, air- 
craft control compensation can be programmed to occur automatically 
when a failure signal is received at the fliglit control comjxiter. This 
control integration feature was assumed in all flights. 

Flight co n tro l compen sa tion effects . - Normal pilot recognition 
delays and reaction to engine failures are of significant concern in 
powered lift aircraft. Cueing by visual means to alert the pilot to en- 
gine thrust loss is one |)otential solution. An alternate approach is for 
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automatic thrust loss comiKMisation pruvidt'd by a thrust command or 
rctmlatcd Uirust control (3). Thrust loss comiKMisation controls such 
as those wore evaluated in the flii'ht simulation program. 

Thrust command '?ontrol: A tlirust control concept is shown in fi|;- 
ure 17. In this concept total eni;ine thrust is compared to the comnuinded 
thrust. After an enitine failure the error nenerated would cause the 
thrust control to advance the |Miwer inputs until the demanded level is 
achievc'd. F^or normal o|K'raiion the pilot can select thrust level sepa- 
rately or topether for all engines. 

Typical flights: Heprcsentative time histones of two flights without 
and with the Uirust command control during an engine failure are shown 
ill figure 18 (withouil and figure 19 (»vith). The loss of the right outlxiard 
engine is simulated during an approach maneuver with subsequent land- 
ing and reverse. Wind is ahead at 15 knots with 0.91 ni|)s (3.0 fps) rnis 
turbulence. Each channel in figures IK and 19 is multiplexed to show 
two variables. Except for the aircraft status data, unfailtKl engine data 
is contained in the ’’short side” traces while tailed engine data is in the 
’’long side” traces. 

Performance without comixMisation: Turbine inlet temiXM’ature 
variations, figure 18, due to pilot throttle movement during the unfailed 
engine approach segment, are small. They were comparable to other 
similar approach flights without engine failure. Variations to III^* K 
(200” H) without substantial ix'aking were typical. 

The outboard engine failure occurs at 263 meters (863 feel) altitude. 
There is a 1.25 second delay in pilot resixinse to initiate a power coni- 
IX'iisation adjustment for the power loss and 4.0 seconds to accomplish 
the adjustment. Roll trim was automatically inserted into the lateral- 
directional stability and command augmentation system at the |xiint of 
failure. As shown on Uie runway jxisition history, the aircraft yaws 
left on tlie approach strip. 

Perfo rma nce with c onnx Misation: Unfailed o|X?ration with the Uirust 
command control in tJie throttle Uuip, figure 19. shows that throttle 
movements to ad|ust thrust on approach tends to produce ix'aking in the 
fuel flow and turbine inlet tenqx'rature. This jx>;Ucing is not observed 
to the same degree in unconi|XMisatc‘d flights. 
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The outboard engine failure occurs at approximately 100 meters 
(325 feet) altitude. Because of the disparity In failure altitudes It Is 
not |X)8sible to draw conclusions regarding thrust command control 
effects on touchdown and landing ix?rformance. At engine failure there 
Is no delay due to pilot reaction. The rapid automatic compensation for 
thrust loss elimuiates the need for pilot corrective action on thrust and 
improves tiirottle resixmse. The controller corrected the power Inputs 
In 3.0 seconds which is about 57 |)ercent of the total time In the un- 
comiiensated case. Roll trim was again automatically engaged. Although 
no judgement can be made In the basis of a single flight, In this case Uie 
aircraft yaw was reduced substantially from the previous flight. The 
|)erformance of Uie unfailed engine after the failure Is very similar to 
normal uncom|x*nsated flights including tlie reverse sequence. 

Witli Uie automatic Uirust command control, peaking occurs In the 
turbine inlet temiierature which n*ay be detrimental to engine life. This 
{leaking is due to the thrust command control which effectively modifies 
the {K)wer demand dynamically even without an engine failure. Off-line 
runs with Uie command system simulated verify Uiis point. The thrust 
command control is not optimized for the system. However, it Is obvi- 
ously interacting with Uie engine control. This points up he require- 
ment for an integrated aircraft propulsion system approach particularly 
for high res{K)nse, close coupled {xiwered lift systems. 

CONCLUDING REMAItKS 

Flight simulation provides a safer and relatively inex{3ensive ap- 
proach to {iropulsion system evaluation coniiiared to flight testing. 

This is es|iecially true for failure modes and effects analysis. And, 
de|)ending on Uie level of detail, it permits observation of internal sys- 
tem {Jaranieters which could oUierwise be difficult or lni)x)sslble to 
measure in flight. 

The real time digital pro{)ulsion system simulation effort under 
the QCSEE Program has generated valuable simulation technology and 
flight simulator operational ex{)erience. It has yielded a technically 
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foaHiblt' rt'ul turn' di(;ital simulation approach which is faster and has 
much higher roal-to-oxccution time ratios than previous nuKlels. The 
techniques develoiH*d have Ihhmi efleciive in iM'rmittini; an excellent level 
of detail within real time constraints. In addition, the application ol this 
nuxlel to a fliiiht simulator lias indicated the need for an intei^ratt'd 
aircraft*pro|Kilsion system approach to |K)wered hit systems. 

Interest in real time di(;ital simulation has grown and apiH'ars to 
have widespread application lu'yoiid flight simulators. These include 
Uie use of the present QCSKK nuKlel in various controls analyses and 
application ol model teciimques develo|M'd in this study to other controls 
analysis models. 
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Figure 3. - Variable pitch fan representation and flow diagram. 
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